Charge fluctuations and electron-phonon interaction in the finite-?/ Hubbard model 
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In this paper we employ a gaussian expansion within the finite-17 slave-bosons formalism to 
investigate the momentum structure of the electron-phonon vertex function in the Hubbard model 
as function of U and n. The suppression of large momentum scattering and the onset a small-q peak 
structure, parametrized by a cut-off Qc, are shown to be essentially ruled by the band narrowing 
factor Zmf due to the electronic correlation. A phase diagram of Zmf and qc in the whole U-n space 
is presented. Our results are in more than qualitative agreement with a recent numerical analysis 
and permit to understand some anomalous features of the Quantum Monte Carlo data. 
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In the past months a variety of experiments have 
pointed out an important role of the electron-phonon 
(el-ph) interaction in many physical properties of the 
cupratesJaSi^ These recent findings have triggered a re- 
newed interest for a theoretical understanding of the el- 
ph properties in strongly correlated systems. 

One of the most remarkable effects of the strong elec- 
tronic correlation on the el-ph properties is a to favor for- 
ward (small q) scattering in the electron-phonon vertex, 
q being the exchanged phonon momentum. This feature 
was investigated in the past by means of analytical tech- 
niques based on slave-bosons or X-operators^i^iSii The 
assumption of forward scattering predominance within 
an el-ph framework was shown to explain in a natu- 
ral way several anomalous properties of cuprates as the 
difference between transport and superconducting el-ph 
coupling constants^ the linear temperature behaviour of 
the resistivity;^ the d-wave symmetry of the supercon- 
ducting gapiiS Small q scattering selection was shown 
moreover to be responsible in a natural way for high-Tc 
superconductivity within the context of the nonadiabatic 
superconductivityiiSi 

So far, this important feature was analyzed only by 
means of analytical approaches in the U — oo limit and 
a definitive confirm of it based on numerical methods was 
lacking. The charge response at finite U was addressed 
in Refs. ITllIT^ which however focused only on q = sus- 
ceptibilities. With these motivations in a recent paper 
Huang et al have addressed this issue in the twodimen- 
sional Hubbard model with generic U by using Quantum 
Monte Carlo (QMC) techniques on a 8 x 8 cluster Their 
results provide a good agreement with the previous an- 
alytical studies and represent an important contribution 
to assess the relevance of el-ph interaction in correlated 
system. 

In this paper we employ the slave-boson techniques 
to investigate the evaluation of the momentum structure 
of the el-ph vertex interaction as function of the Hub- 
bard repulsion U . While the previous analytical studies 
were limited to the U = oo limit ^i^iSii we were able in 
this way to evaluate the small momenta selection in the 



whole phase diagram of parameters U-n, where n is the 
electron filling. In particular we show that the predomi- 
nance of small q scattering is strongly dependent on the 
closeness of the system to the metal-insulator transition 
(MIT). Our analysis suggests that forward scattering is 
essentially ruled by the one-particle (mean-field) param- 
eter Zmf which represents the band narrowing factor due 
to the correlation effects. We also discuss some anoma- 
lous features of the numerical results of Ref. [I^- The 
fair agreement between our slave-boson results and the 
Quantum Monte Carlo datai^ supports the reliability of 
our approach to investigate these issues. 

Working tool of this paper will be the four slave-boson 
method first introduces by Kotliar and Ruckenstein to 
investigate the Hubbard model at finite U. In this ap- 
proach the physical Hilbert space is enlarged by intro- 
ducing, in addition to the standard fermionic degrees of 
freedom Cicr, four auxiliary bosons, e, d Pa-, which count 
respectively empty, double and single occupied (with spin 
a) sites The extended Hilbert space is thus restricted 

to the physical one by the Lagrange multipliers A^^^ , Aa^'' . 
The Hubbard Hamiltonian reads thus: 
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where tij is the bare hopping term, here considered 
within a nearest neighbor model, and Z[. . .] is the reduc- 
tion of the hopping term due to the correlation effects. 
A standard method to account for charge, spin, and dif- 
ferent kinds of fluctuations is to expand Eq. ^ at the 
gaussian level around its mean-field solution for the aux- 
ihary fields e, d, p^, A^^), A^^^iiJ^ Eq. Q can be thus 
re- written as: 
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where the first term describes the mean-field approxima- 
tion where all the auxiliary fields are determined by their 
self-consistent saddle point solution, the second term are 
the quadratic fluctuations of the auxiliary fields about 
their mean-field solution, comprehensive of the electronic 
loop self-energy contributionSfii*i^ and the third one is 
the linear electron-slave-boson interaction. The index /i, 
V run on the seven auxiliary fields. 

The el-ph vertex function g{k.^ q, uj) in the presence of 
electronic correlation can be promptly determined as the 
linear response of the systems described by Eq. ^ to 
an external charge excitation. In the paramagnetic case 
the spin index a can be dropped and the auxiliary field 
space reduce to five components. At the gaussian level 
we haveiiii^ 

5(k,q;..) = 1 -2Y,K,^ [D-^u)]"" (3) 
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where the factor 2 takes into account the spin degen- 
eracy, is the renormalized slave-boson propagator, 
2 [D-Hw)]^" = 2 [Sqi]"" + n^''(tj), which includes the 
"bubble" contribution to the slave-boson self-energy 

/(i?k) - /(i^k+q) 
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and where 
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Here f{x) is the Fermi factor and i?k is the renormalized 
electronic dispersion Ey^ = ^MF^k- 

We are now in the position to apply Eqs. ©-© to cal- 
culate the q-structure of the electron- phonon vertex func- 
tion in the presence of electronic correlation for generic 
Hubbard repulsion U and electron filling n. From the 
technical point of view we use a slightly different ver- 
sion of the four slave-boson technique as introduced by 
Kotliar and Ruckentstein. We note indeed that Eq. ||2J), 
when derived in the context of the Kotliar-Ruckenstein 
formalism, present serious pathologies at finite U in the 
limit n — > 0, 2 where the linear electron-slave-boson ma- 
trix elements diverge«i& As pointed out in Ref. how- 
ever several choices for the auxiliary fields e, d, pi^j can be 
employed as long as they fulfill the constraints imposed 
by the real Hilbert space. We make use of this freedom 
to avoid those unphysical divergences. The analytical 
derivation of Eq. ^ in this new basis of auxiliary fields 
is long and cumbersome, but quite straightforward and 
it does not affect in a significant way our results which 
focus on the MIT transition Technical details will be thus 
presented in a longer publication. 

In Fig. plot the evolution of momentum structure 
of the electron phonon vertex function g{k-p^ q; w = 0) for 



FIG. 1: (color online) Momentum dependence of the electron- 
phonon vertex function g{kF,c\;,uj — 0) in the presence of 
electronic correlation for different U/Uc and n. In each 
panel solid lines represent, from the top to the bottom: 
n = 0.1, 0.3, 0.5, 0.7, 0.9. Standard labels for the high symme- 
try points correspond to F = (0,0), X = (0,7r), M — {iT,n). 



UJ = and k = /cp lying on the Fermi surface for different 
electron filling n and U. The top panel in Fig. ^ repro- 
duces the n-dependence of ^(fcF, q; a; = 0) in the large-?7 
limit. It agrees quite well with the previous studies of 
Refs. 00 and it shows the the suppression of the large- 
q scattering as the MIT is approaches for n ^ 1 and 
the corresponding development of a small-q momentum 
structure. As a marginal difference we can note that 
the four-slave-boson approach at finite U predicts a less 
pronounced peaked structure than the U = oo one slave- 
boson method. This slight discrepancy can be traced 
back to the different band narrowing factor between our 
case {Zmf ~ 2(5/(1 -I- S) as in Ref. [ij and the one-slave- 
boson U = oo treatment (Zmf ^ <5). 

One of the main results of our analysis is that it per- 
mits to investigate the MIT and the momentum structure 
of the el-ph vertex along a different path, namely increas- 
ing U for fixed n. As shown by the comparison between 
the three panels of Fig. ^ the predominance of forward 
scattering is effectively suppressed as long as we move far 
from the MIT by reducing the Hubbard repulsion J7, in 
similar way with what is found as function of n. This 
observation suggests that the small-q structure is only 
ruled by some macroscopic parameter which quantifies 
the closeness to the metal-insulator transition, regardless 
the microscopic parameters as U or n. A natural candi- 
date along this perspective is the band narrowing Zmf, 
which vanishes Zmf — > at the MIT (n = 1, ?7 > Uc) and 
Zmf ^ 1 in the uncorrelated limit (n = 0, 2 or [/ = 0). 

The effectiveness of the electronic correlation in giving 
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FIG. 2: (color online) Plot of the isolines for the quantities 
Zmf (panel a) and Qc (panel b) in the U-n space. Since 
both quantities are symmetric for reflection with respect to 
the n = 2 axis, they are shown respectively in just one 
half-plane. Solid lines are (panel a, from right to the left): 
Zf = 0.1,0.2,0.4,0.6,0.8,0.9; (panel b, from left to right): 
Qc/tt = 0.2, 0.3, 0.4, 0.5, 0.6, 0.8, 1.0, 1.2, 1.4. Panel (c): plot of 
Qc VS. Zmf for different electron filling n. 



rise to a small-q structure in the electron-phonon scatter- 
ing has been previously discussed in literature in terms 
of a cut-off parameter On the physical ground qc is 
expected to be related to the inverse of the correlation 
length ^, qc oc 1/^. At the MIT ^ ^ oo and 0. 
For the practical purposes, it is often difficult to extract 
a cut-off parameter qc from a complex structure as those 
shown in Fig. ^ To this aim we define a cut-off qc as the 
full- width at the half-maximum (FWHM) of the electron- 
phonon vertex g(k-p,ci;uj = 0) along the T-M line with 
respect its q = value. We would like to stress that this 
definition of qc is just one of the possible choices and it 
does not pretend to be rigorous. As matter of facts for 
systems far from the MIT, where no well defined q = 
peak is recovered, no qc according such a definition can 
be found. This should be interpreted as the absence of 
significant small-q modulations induced by the electronic 
correlation. Note, on the other hand, that close to the 
MIT the small-q peak can be always approximated as a 
Lorentzian, and the definition of qc introduced above is 
indeed directly related to the inverse of the correlation 
length ^. 

In Fig. I^t-jb we plot the Zmf- and gc-isolines as ob- 



tained by applying our FWHM definition in the whole 
U-n phase diagram. The striking similarities between 
panels (a) and (b) confirm that the small-q structure is 
mainly ruled by the strong band narrowing close to the 
MIT. This is pointed out in the most remarkable way 
in panel (c) where the cut-off qc, evaluated for different 
n and different U, is shown to lie on an almost univer- 
sal curve when plotted as function of Zmf- This is a 
quite important result because the band narrowing fac- 
tor Zmf, as well as quasi-particle spectral weight Zqp 
which is strongly related to Zmf, is in principle a physi- 
cal quantity which can be experimentally determined to 
provide an estimate of qc- The universality of qc as func- 
tion of Zmf points also towards some general mechanism 
responsible for the small-q scattering. This mechanism 
has been discussed in literature in terms of proximity to 
a phase separation which reflects in a charge response in- 
stability at q = 0. Before addressing the connection be- 
tween small-q scattering and phase separation, we would 
note however that, although the universal behaviour of 
qc vs. Zmf appears quite robust, some discrepancies ap- 
pear for values of n very close to n — 1. This is also 
in agreement with a comparison between panels (a) and 
(b) which shows that for n = 1, in contrast to the Zmf 
isolines, all the gc-curves converge to the U ~ Uc point, 
where qc jumps from qc ~ t: to qc ~ from U — U~ to 
U — U^. A closer look at this discrepancy shows that 
it can be related to the presence of the Van Hove sin- 
gularity at n = 1. Indeed, while the quantity Zmf does 
not depend on the details of the electronic dispersion, 
in this case the q = response functions are described 
by an ideal metal with infinite density of states for any 
U < Uc ■ It can be checked that the introduction of a 
finite t' removes this anomalous features and qc evolves 
in a smooth way increasing U from [/ = to f7 = J7c for 
n = 1. 

As a final issue of this Letter we would like to com- 
ment about the robustness of our results, based on a 
gaussian expansion within the context of a finite-?/ slave- 
boson formalism, with respect the inclusion of higher or- 
der terms. A good check in this perspective would be the 
comparison of our data with q-resolved numerical calcu- 
lations. Unfortunately, simulation and numerical work 
along this line is lacking in literature since the widely dif- 
fuse dynamical mean-field approach mainly probes only 
local (-integrated) quantities whereas exact diagonaliza- 
tion and QMC techniques are limited by finite size of the 
cluster and by finite temperature effects. A seminal step 
to fill the gap has been provided by Huang et al. who 
provided in Ref. .13i . for the first time at our knowledge, 
numerical support by QMC simulations for the insur- 
gence of small-q structure in the el-ph vertex function by 
increasing U . As a by-product of their analysis, Huang et 
al. reported also an interesting increasing of the absolute 
magnitude of the el-ph vertex function for small momen- 
tum q = (7r/4, 7r/4) in the large U regime, whereas the 
increase is absent for large momentum q = (tt, tt). This 
is quite puzzling because the development of a small-q 
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FIG. 3: (color online) (a) Dependence of the el-ph vertex 
function on U for n = 0.88 and for q = (7r/4, 7r/4) (solid line, 
/3 = 2), q = (tt.tt) (dashed line, f3 = 2), and q = (7r/4, 7r/4) 
(dot-dashed line, /3 = 50). (b) Phase diagram for phase sep- 
aration (q = 0) and charge instability (q = 7r/4, 7r/4). The 
dashed line shows the factor Zmf evaluated along the phase 
separation instability line 



structure is thought to be relative to the large momen- 
tum scattering and it is usually not accompanied with 
an increase of the magnitude of g at q = 0. The origin 
of this anomalous feature has been not well understood 
so far and a possible connection with charge excitations 
driven by the exchange term J oc t'^/U in the large U 
limit was mentioned 

In order to test to reliability of our results we have ad- 
dressed this issue within the context of our slave-boson 
formalism. In Fig. we show the dependence of the el- 
ph vertex function g{kF, q; a; = 0) on the Hubbard repul- 
sion U as calculated within our approach for the same pa- 
rameters as in Ref. 13, namely n — 0.88, q — {tt / 4, n / A) , 
q = (tt, tt), /3 = 2, where /? is the inverse of the tem- 
perature T in units of the nearest-neighbor hopping el- 
ement t. The agreement between our findings and the 
data of Ref. is quite excellent, although quantitative 
differences are of course present. In particular we find 
that, while for large momenta f/(fcF,q;u; = 0) steadily 
decreases with U, for q — {tt/A, 7r/4) it shows an upturn 



for J7 ~ 8 until a divergence occurs for Uc — 9.3, sig- 
nalizing an charge instabilityi-i According this view one 
is attempted to associate the upturn of the el-ph ver- 
tex function g{kp, q;uj = 0) as function of U reported in 
Ref. ^3 3.S an incipient transition towards some charge 
instability or phase separation^ This does not imply 
however that phase separation is effectively established, 
and it should be remarked that the actual occurrence of 
phase separation in the Hubbard model is still object of 
debateii^ In particular the large temperature /3 = 2 at 
which the QMC were performed could play an crucial 
role in that. To investigate this final issue we evaluated 
again the el-ph vertex function for the same parameters 
but a much lower temperature /3 = 50. As shown in Fig. 

quite surprising the charge instability disappears de- 
creasing the temperature. 



In order to understand in more detail the origin of the 
phase separation instability as function of the tempera- 
ture T we show in Fig. |3|d the phase diagram in the T 
vs. U space with respect to phase separation (q = 0) and 
to charge ordering (q = (7r/4, 7r/4)). Phase separation 
as well as charge instabilities occur in our slave-boson 
calculations only above a certain temperature T/t ^ 0.2 
(/3 ^ 5). A interesting insight comes from the compari- 
son of the critical temperature Tc at which the instability 
towards phase separation occurs with the band narrow- 
ing factor Zmf (Fig- El- The similar dependence of Tc 
and Zmf on U points out that the onset of phase separa- 
tion by increasing temperature is ruled by the comparison 
between the temperature T and the effective bandwidth 
W — ZuF^t energy scales. In particular phase separa- 
tion is established when T becomes of the same order of 
ZuFt- We would like to stress that in principle the phase 
separation instability found by our slave-boson calcula- 
tions could be washed out when higher order fluctuations 
than the gaussian ones are taken into account. The agree- 
ment between our data and QMC results suggests how- 
ever that the correlation between small-q modulation, 
phase separation and the upturn as function of U could 
be a robust feature of correlated systems. At a specula- 
tive level, we think that the small coherent part of the 
electronic spectral weight (with the electronic dispersion 
scaling with Zmf) could be responsible for the tendency 
towards phase separation, which is however prevent by 
the presence of the incoherent states. Further analytical 
and numerical work on this direction will help to shed 
new light on this issue. 
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